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Abstract 

In this paper we focus on the large n probability distribution function of the largest 
eigenvalue in the Gaussian Orthogonal Ensemble ofnxn matrices (GOE n ). We prove 
an Edgeworth type Theorem for the largest eigenvalue probability distribution function of 
GOE n . The correction terms to the limiting probability distribution are expressed in terms 
of the same Painleve II functions appearing in the Tracy- Widom distribution. We conclude 
with a brief discussion of the GSE n case. 



1 Introduction 

Limiting probability distributions laws from Random Matrix Theory have found many 
applications outside their initial domain of discovery; the length of the longest in- 
creasing subsequence ( P. Deift et al. [1]) properly scaled converges in distribution to 
the Unitary Tracy- Widom law, the properly scaled largest principal component of a 
white Wishart converges in distribution to the Orthogonal Tracy- Widom law (I. M. 
Johnstone [H] ) . For recent reviews we refer the reader to [U [5j El [131 12Z] • In these 
applications of Tracy- Widom distributions, one would like to control quantitatively 
the range of validity of the various limit laws. One therefore needs finite n correction 
to these limiting distributions. In a previous work [2] we initiated the study of this 
problem for the Gaussian Unitary Ensemble of n by n matrices (GUE„). Following 
this work, we will derive the analogous result for the Gaussian Orthogonal Ensemble 
(GOE n ) in this paper. The derivation of the probability distribution function of the 
largest eigenvalue for Gaussian Symplectic Ensemble (GSE n ) is similar to the GOE ?t 
case up to the parity of the size n of matrices in consideration. We will therefore 
mention at each step of the derivation the corresponding result without much expla- 
nations except when it is necessary. We seek a large n expansion of the probability 
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distribution function of the largest eigenvalue in GOE n and GSE n similar to the 
following Edgeworth Expansion arising in probability in applications of the Central 
Limit Theorem. 

If S n is a sum of i.i.d. random variables Xj, each with mean /i and variance a 2 , 
then the distribution F n of the normalized random variable {S n — nii)/{a^/n) satisfies 
the Edgeworth expansion^ 

r 

F n (x)-$(x) = (j)(x)J2 n ~^ +lR A x ) + o(n-^ r+1 ) (1.1) 

i=3 

uniformly in x. Here $ is the standard normal distribution with density <ft, and Rj 
are polynomials depending only on E(X*) but not on n and r (or the underlying 
distribution of the Xj ) . 

If we view the random matrix ensembles of n by n matrices in terms of the 
associated eigenvalues, then the Gaussian /3-ensembles are probability spaces on n- 
tuples of random variables {Ai, . . . , A n } (think of them as eigenvalues of a randomly 
chosen matrix from the ensemble.) The probability density that the eigenvalues lie 
in an infinitesimal intervals about the points x\, . . . ,x n is 

P n/3 (xi, • • • ,x n ) = C raj9 exp ( ~ ) TT \ x i ~ x kf, (I- 2 ) 

V i / j<k 

with 

— oo < Xi < oo, for i = 1, • • • , n. (1.3) 

Here C n p is the normalizing constant such that the total integral over the x^s is one. 
The cases (3 — 1, 2, 4 correspond to the GOE n , GUE„ and GSE n respectively. We 
denote the largest eigenvalue by \ Max , and 

F n>p {t) = P(A^ < t) (1.4) 

the probability distribution function. 

When (3 = 2, the harmonic oscillator wave functions 

, , 1 



H k {x) e~ x /2 k = 0, 1, 2, . . . 



obtained by orthonormalizing the sequence x k e x2 (with H^ix) the Hermite polyno- 
mials of degree k) play an important role. We also have the Hermite kernel 

v I \ / \ / \ <Pn{x)(p n -l(y) ~ <fn(y)<fn-l(x) 

K ni2 {x : y) = 2_^<p k (x)cp k {y) = J- , (1.5) 

fc=o V 1 x y 

which is the kernel of an integral operator K 2 acting on L 2 (t, oo) , with resolvent 

R n (x,y-t) = (I — K n y 1 K n (x, y). (1.6) 



1 We assume, of course, the moments E(Xj c ), k = 3,...,r, exist; and as well, the condition 
lim^i^oo sup \ip(Q\ < oo where ip is the characteristic function of Xj, see [7]. 



2 



The product on the right is operator multiplication. We have the following represen- 
tation of (11.21) . (see for example, 



P n2 0i, • • • ,x n ) = det(i^ ri (x i ,x j ))i< ij < n . 
Following Tracy and Widom [21], we define 



/ n \ 4 / n \ 4 

= ( 2 J <Pn(x), 1p(x) = I - J tp n -l(x), (1.7) 

by e the integral operator with kernel 

e t (x) = -sgn(x-t), (1.8) 
P the differentiation with respect to the independent variable, 

Q n>i (x;t) = ( (/ - Kn)- 1 , xVn) (1-9) 

and 

P ft>i (z; *) = ((/- An) -1 , xVn-i). (1.10) 
We introduce the following quantities 

= Qn,ifct), Pn ,i(t) = PndW) ( L11 ) 

«n,i(*) = {Qn,h<Pn), V n ^(t) = {P n ,i,<fn), (1-12) 

«n,i(<) = (Qn,i, Vn-l), and W n ,i(<) = (■Pn.tj^n-l)- (1-13) 

Here (•,■) denotes the inner product on L 2 (t, oo). In our notation, the subscript 
without the n represents the scaled limit of that quantity when n goes to infinity, and 
we dropped the second subscript i when it is zero. 

If Ai is Airy function, the kernel K n> <i(x, y) then scales^] to the Airy kernel 

K M (X,Y) = A-WAi'OT-AiOTAi'W (L14) 

Our conventions follows: 

Q l (x;s) = ((I -Km)' 1 , x*Ai), Q (x;s) = Q(x;s), (1.15) 

P(x; s) = ( (I — A^Ai)^ 1 , x l Ai'), P (x; s) = P(x; s), (1.16) 

?i(s) = go(s) = g(s), p t (s) = P(s;s), Po(*)=P(s), (1.17) 

«i(s) = (Qi, Ai), u (s) = u(s), Vi(s) = (P i5 Ai), u (s) = u(s), (1.18) 
^i( s ) = (Qi, Ai'), v (s)=v(s), Wi(t) = (Pi, Ai'), and w {s)=w(s). (1.19) 



2 For the precise definition of this scaling, see the next section 
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Here ( • , • ) denotes the inner product on L 2 (s, oo) and % = 0, 1, 2. 

We also note that q(s) is the solution to the Paileve II equation q"(s) = sq(s) +2g 3 (s) 

with the boundary condition q(s) ~ Ai(s) as s — > oo. 

We use the subscript n for unsealed quantities only. 

ll n>1 := I Rn(x,t;t)dx, V n ,i ■= / P n (x;t)dx, Q n ,i ■= / Q n (x;t)dx, 

J — oo J — oo J — oo 

(1.20) 

and 

/OO /'DC 
e t (x)Rn(x,t;i)dx, V nA (t) := / e 4 (x)P n (x; t)dx, 
-oo J —oo 

/oo 

£ i (x)g n (x;t)da;. (1.21) 
■oo 

The epsilon quantities are 

Q„, e (x;t) = ( (/ — K n )~ x {x, y) , Sip(yj), q n , e {t) = Q n , e (t;t) (1.22) 

VW = {QnA^t)^^)), «n,e(*) = (<5n, £ (^ *)> ^(^)) , (!- 2 3) 

where ( ■ , ■ ) denotes the inner product on L 2 (t, oo). 

The GOE n and GSE n analogue of (11.281) in Theorem 11.11 bellow will follow from 
representations (equations (40) and (41) of [24"].) 



Fn,i(t) 2 = F n>2 (t) ■ ((l - v n>E {t)) (1 - ^ n>1 (t)) - ^(q n>£ {t) - c^Vn^tYj (1.24) 



and 



F nA {t/V2f = F n>2 {t) -((I- v n>£ {t)) (1 + ^n nA (t)) + ^q n>£ {t) V nA (t^j . (1.25) 

Here we first derive a large n-expansion of TZ n ,i, V n ,i, 7^n,4> V n ^ v n ,e, andg nj£ in terms 
of p n and q n , by solving the associated systems of differential equations. We then 
substitute the resulting expressions in (11.24ft and (11.251) . We will need the following 
result which gives the large n expansion of F n , 2 . 

Theorem 1.1. If we set 

t = (2(n + c))3 + 2~^n~^ s and (1.26) 
E C:2 := E Cj2 (s) = 2wi - 3u 2 + (-20c 2 + 3)v + uiv® - u vi + u Vq - u&dq. (1.27) 
Then as n — > oo 

F n>2 (t) = F 2 (s) |l + cu (s)n-^ - l^ Ci2 ( s )n-i| + 0{n~ x ) (1.28) 

uniformly in s, and 

F 2 (s) = hmF nj2 (t) = exp(- f (x - s)q(x) 2 dx] (1.29) 
n ^°° V Js J 

is the Tracy-Widom distribution. 
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1.1 Statement of our results 

To state our main result we need the following definitions 



a := a (s 



q(x) u(x) dx, 

/oo 
q(x)dx, 



v := v{s) 



p(x)dx = a(s) — q(s), 



1.30) 

1.31) 
1.32) 



7] := 7](s) 



1 f°° 

—= / (6qv + 3pu + 2p 2 + 2piv + 2pv\ — 2q 2 u — 2qiU\ — 2qu 2 ) (x) dx 
20v2 Js 

20c 2 q'{s) + 3p(s) 



E c , 2 (s) e 



20^ 



ca cp (2c — 1) i/ 2 



+ 



20 ^~ J ~ 2/x 2 2/i 
, u(u + 8cq) 7] 



v 



+ c«^--(l-e-": 



1.33) 



+ 



32,U 



4v/2 



cV 2 - u 



( 2V2cV-3r/ v 2 -8{2cp + c 2 q 2 ) - Ac 2 a 2 
+ 32^ 

-(4 C V + 3cV-^ 2 )^^. (1-34) 
v 7 8p 2 



Unfortunately, we did not find a simple representation of t] and E c \. Nevertheless 
the quantities a(s), fi(s), v(s), rj(s) and E c> i(s) are easy to compute numerically. For 
E C) i(s) and rj(s) we only need the recurrences relations defining Pi(s) and qi(s) in 
term of q(s) and q'(s). We find the following representation of the large n probability 
distribution function for X max in GOE n . The derivation of the analogous result for 
the GSE n follows exactly the same steps as the one given in this paper. Here is our 
main result. 



Theorem 1.2. We set 



Then as n — > oo 
F n ,i{tf = F 2 (s) 



t = (2(n + c))a +2 2fl as. 



1.35) 



-fj,(s) 



iq(s) + u(s))e-^-^(l-e-^) 



n a + 



1.36) 



uniformly for bounded s. 
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Note that unlike the GUE n case where for c = the n~3 correction term vanishes 
as shown in equation (ll.28p . the correction term does not vanish in the GOE„ 
no matter what the fine tuning constant c is. 

In §2 we reproduce the derivation of (11.241) following Tracy and Widom in [23]. In 
§3 we solve the system of equations satisfied by the various functions on the right of 
(ll.24p for our derivation of the Edgeworth expansion of the probability distribution 
of the largest eigenvalue in the GOE n . 



2 Derivation of F n p 

We treat here the case n even. If we set 

K . ( K n! 2 + ip®e(p K n! 2 D - ip ® ip \ . . 

M ' V eK n>2 -e + eij)®etp K n>2 + e^p ® <p J { ' } 

then F% is the Fredholm determinant of K n l on the set J = (t, oo). If we denote by 
X the multiplication by the function Xj( x )j then F^^t) is the Fredholm determinant 
of the integral operator with kerneH 

K i = yAx) ( K n,2 + ^®e V K n , 2 D -ij®ip\ , 
nA XJ[ ] \ eK n>2 -e + eiP®ecp K n>2 + eiP®<p J XAV} 

( K n,2 + ^®SLp K ns2 D — Ip ® ip \ . . 

" X \ e K n>2 -e + eip®etp K n>2 + etp ® <p ) X 1 ' ' 

on M , see for example [23] equation (31). Using the following commutators, 

[K n>2 , D}=tp®il) + il)®tp, [e, K n p] = -eip <g> eip - eip <g> e<p (2.3) 

(ip and (f appear as a consequence of the Christophel Darboux formula applied to 
K n> 2,) we have 

K n ,2 + ip ® <p = D sK n! 2 + Deip ®eip = D(eK nt2 + eip ® eip) = D(K nfl e - ep <g> eip) 
K n ,2 D - ip <g> e = D K n>2 + p®ip = DK 2 + Deip <g> e = D (K n>2 + eip ® -0) 

and 

^ ^n,2 — £ + £ ^ ® ^V 9 = ^n,2 £ ~ £ ~ e (p ® eip 

as .De = /. Our kernel is now 

K v ( D ( Kn ' 2 £ - £< f®^) D (K n>2 + e<p®iP)\ . . 

x I K ni2 e-e-eip®eip K n , 2 + eiP®p J x { ' 



(2.5) 



3 To simplify notation we kept the same notation for the integral operator as well as the kernel 
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Since K n \ is of the form AB, we can use the fact that det(J — AB) = det (7 — BA) 
and deduce that the Fredholm determinant of K n \ is unchanged if instead we take 
K n i to be 

(K nj2 e-etp®eip)x (K n , 2 + etp ® ip)x \ ( X D \ . . 
(K nj2 e-e-e<p®eip)x (K nft + <g> <p) x J ' \ X J 

(K n , 2 e-e(p®eip)xD (K n , 2 + etp ® ip)x \ ? 2 7 ^ 

(K n , 2 e-e-eip®etp)xD (K U:2 + etp <g> </>) x / 



det(7 - i^ n ,i) = det 



/ - (iT n;2 £ - Elf ® £ljj)x D -(#n,2 + ^ ® VOX 

-(K n ,2 s-s-eip®eip)xD I - {K n>2 + eip<g>tp)x 

n (2 ' 8) 

Performing row and column operations on the matrbQ does not change the Fredholm 
determinant. We first subtract row 1 from row 2, next we add column 2 to column 1 
to have the following matrix 



I - {Kn t2 E-etp® eip)xD - (K nt2 + etp® ip)x ~{K n>2 + etp® ip)x 
e X D I 



(2.9) 



Right- multiply column 2 by —exD and add it to column 1, and multiply row 2 by 

(Kn,2 + ® i>)x an d a dd it to row 1 to have 

I - (K nf2 e- etp®Eip)xD- {K Ui2 + etp ® ip)x + (K nt2 + e<p ® ip)x£X D 

I 

(2.10) 

We therefore have, 

det (J - K nA ) = det(l -(K ni2 £-e<p®£i/j)xD + (K n , 2 + e<p®<ip)x{ExD-l)) 

(2.11) 

= det (/ - K n . 2X - K n , 2 (I - x)exD- (e<p ® V)(x - - ^ ® ^e X D) (2.12) 

We used the fact that e is antisymmetric to have 

E(f ® Elp X D = E(f ® tpE 1 x D — —sip ®tjjExD, 

and if we note that x is multiplication, then the determinant is 

= det (I - K n>2 x - K n>2 (I - x)£XD - eip®ip(l- x)^xD - ep ® x$) ■ 

Next we factor out / — K n 2 and note that (7 — K n ^)~ x = I + R n ,2, where R n ,2 was 
defined as the resolvent of K n>2 , and (I — K n ^)~ x Eip = Q n ,e- We are interested on the 
determinant of the following operator 

(l-K n!2X )(l-(K n>2 +R n , 2 K n>2 )(I-x)£xD-Q n ,^ (2.13) 

We have a large n-expansion of det (I — K Uj2 x) = Fn.2 from the author work in 
GUE n see [2] or [3]. We will therefore focus our attention on the second factor of 



4 This does not change the determinant, for more details see [24] 
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f)2.13p . We will represent this factor in the form (/ — Y^j=\ a j ® Pj) and use the 
well known formula det(7 — X^=i a j ® Pj) — det(5jj — (atj, Pj)) i j =1 k to expand the 
Fredholm determinant. First we need to find a representation of ex D as a finite rank 
operator. To this end we introduce in this section the following notation, 

e k {x) = e(x-a k ), R k (x) = R ni2 (x,a k ), 8 k (x) = 8(x-a k ), ai=t, and a 2 = oo. 

With the new notation J = (t, oo) = (ai, 02), and the commutator 

[x £>] = -£1 ® 5i + 6 2 ® 5 2 , 

gives 

e[x -D] = -£ X <g> 5i + e 2 ® ^2- 
Next we use the identity e D — I to have 

(J - X )£X£> = (J - X )e [X D] = -(I - X )e 1 <g> 5i + (/ - x)e 2 ® fc, (2-14) 
and the representation 
{K n>2 + R n , 2 K n , 2 )(I -x)e X D=J2 (- l T( K n,2 + R n .2K n , 2 ){I - x )e k ® <f fc . (2.15) 

fc=l,2 

We substitute (12.141) and (12. 15ft in the second factor of ( 12 . 131) and have, 

I- (- l ) k ( K n,2+Rn,2K n i)(I-x)e k ®5 k - (-^Qn^ ' {l-X)^k®h-Qn,e®X^ ■ 

k=l.2 k=l.2 

(2.16) 

The dot in this formula represent operator multiplication. In this case we just multiply 
the kernels using the formula (a (g) /3)(7<8> 5) = (j3, 7)0 (g) 5, to have the following form 

of (EUD 

I- (- l ) k ( K n,2+Rn,2K n)2 )(I-X)e k ®5 k - ^ (^~XH) Qn,e^ k -Qn,e®X^ ■ 

k=l,2 k=l,2 

(2.17) 

We have 

£2 = --, (J-x)ei = (J-x)^2 = ~, and i2 2 = i^afa, 00) = 0. 
If we substitute these value in (12.171) . it then becomes, 

I ~ Qn,e ®X^~\ [(K,2 + Rn,2 K n , 2 )(I - X ) + (V, V ~ X))Qn,e] ® (<*1 " 5 2 ). (2.18) 

This operator is of the desired form 

I - ^2 a k ®(3 k 

k=l,2 
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with 

(Xl = Qn,e, <*2 = 2 [(kn,2 + Rn,2Kn,2)(I-X) + aiQn,e\, Pi = X&> p2 = $l-$2, (2.19) 

and 

o-x = {i'A 1 - x))- 

The corresponding inner product are; 

= v n>£ , (cki, f3 2 ) = q n ,e + c v (2.20) 

with 



oc 



If 1 

c 9 = e <p{oo) = - J ip(x)dx c^p = etp(oo) = - I tp(x)dx, (2.21) 

and for n even 



and 



c, = (^^) 1/4 2- 3/4 - n/2 |^ 2 y, (2-22) 



(a 2 ,/?i) = ^CPn,i - ai + aiu nie ), (a 2 , /3 2 ) = ^(7^ n ,i + ai<?n, e - aic v ). (2.23) 
The determinant of (12.171) is therefore 

(1 - v e )(l - -K n , x ) - -{q n , e - Cip )V nA . (2.24) 

In a similar way we obtain the second factor on the right side of ( I1.25P for the GSE„ 
case 

(1 - v n , E ) (1 + 2^n, 4 ) + -QnAt) V nA {t). (2.25) 

We will derive differential equations involving the various terms in equation (12.241) 
and solutions in terms of q n and p n . Then used known asymptotic of q n and p n for 
the derivation of a large n expansion of F n \. 



3 Differential Equations 

3.1 System of Differential Equations 

This section will be devoted to finding expressions of 

T^n,i(t), V n ,i(t), q n ,e(t), and v ntE (t), 
and show how to obtain the corresponding quantities 

TZ nA (t), and V nA (t) 



9 



for the GSE n . 

To solve the associated system of differential equations it is convenient to introduce 
the following quantities 

and 

Pn,2(x,y) the kernel of (J - K n ^xY x ■ (3.1) 

We have 

Pn,2{x,y) = S( x ~y) + Rn,2{x,y), (3.2) 

and 

d d f l f l d 

—TZ n ,i (*) = 37 / R n .2{x,t)dx = R nt2 (t,t) + —R n2 (x,t)dx. (3.3) 
dt dt ' dt 

Formula (45) of [21] gives 

^-Rn^ix, t) = -—Rnpfa t) - p n (t)Q n (x; t) - q n P n {x; t), (3.4) 
dt dx 



and we find that 



<,i(t) = J t n n,i{t) = -p n (t) Q„,i(t) - q n (t) V n>1 {t). (3.5) 



We also have 

d „ d r* 



and 



Qn,i(*) = ^Sn,i(*) = ^/ Qn(x;t)dx = q n (t)(l-n nA (t)), (3.6) 

= Tt Vn ' 1 ^ = It J Pn ^ dx = Pn(t)(i-n njl (t)), (3.7) 

where we used 

d d 

—Q n (x-,t) = -q n (t)R n2 (x,t), and —P n (x;t) = -p n (t)R n 2(x,t). (3.8) 
dt dt 

The other derivatives are 

d d f°° f°° d 

^ t Un > £ ^ = ~d t J Qn(x,t)eip(x)dx = -q n (t)e<p(t) + J — Q n (x,t)e <p(x) dx (3.9) 

/OO ^ /*oo 
— i^OM)^^)^) = ~9n(*) J Pn,2(%, t)e <p(x) dx, (3.10) 

and therefore 

«n,e(*) = -Qni^qnAt)- ( 3 -H) 

Similarly 

d „ d r 00 f°° d 

~jf>n,e{t) = ^ y P n (x,t)e(p(x)dx = -p n (t)e(p(t)+j —P n (x,t)e(p(x)dx, (3.12) 
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or 

v 



i t£ (t) = -p n {t)q n , £ {t). (3.13) 
The last of these is 

d d f 

y) ey(y)dy - g„(t) (xPn(y; t),etp(y)) - p n (t) (xQ n (y; t),e<p(y)). 



d 

d~y Pn - 



Integration by parts together with the boundary conditions and Ds = I gives 
r P> r 



8_ 

dy 



Pn,2{t,y)e<p(y)dy= J p n>2 (t,y)x£v(y)dy = q n {t), (3.16) 
which in turn gives 

= Qn{t) ~ v n , £ (t)q n (t) - u n , £ (t)q n (t) . (3.17) 
The boundary conditions at t = oo for these function are, 

7^n,i(oo) = 0, <2 n ,i(oo) = 2c v and V n ,i(oo) = 2c^ = as n is even. (3.18) 
and 

Vn,e{°°) = °> u n>£ (oo) = 0, and g„, e (oo) = c v . (3.19) 
The associated systems of equations are; 

?n,e(*) = ?«(*)(!- 5 n,e(*)) - Pn(*) 

(1 -v n ,e)'{t) = Pn (t)q n , £ (t); (3.20) 

«n, 6 (*) = -Qn{t)q n ,e{t)- 



and 



(l-ft n ,i)'(t) = Vn{t) Qn,l{t) + <ln{t)V n ,l{t)\ 

Q' n ,i(t) = q n (t)(l-TZ nA (t)); (3.21) 

Kiit) = Pn (t)(i-n n ^t)). 

3.2 Asymptotic solutions 

We will define in this subsection only 

V n , E {t) = 1 - v n>£ {t), and TZ n>1 (t) = 1 - TZ njl (t). (3.22) 

With this notation, system (13.201) is 

, / u n , £ (t) \ ( -g n (*) \ / u n>£ (t) \ 

-r, V n . £ (t) = p n (t) ■ V n>e (t) (3.23) 
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and fl3T2Tj) is 

Qn,iO) \ / q n (t) \ ( Q n ,i(t) 

V n ,x(t) = Pn (t) ■ V n>1 (t) (3.24) 

n n ,l{t) ) V Pn(t) Qn(t) J \ K n>1 (t) 

If we let 

X* = (u n , E (t),V n>E (t),q n , E (t)) and Y l = (Q n ,i(t), P„,i(t), ^„,i(t)) , 
then (13.201) and (I3.2ip have the following representations 

X\t) = A(t)X(t) and Y' (t) = B(t)Y{t) (3.25) 

with 

X\oo) = (0,1, c^), and F'(oo) = (2c v , 0, 1). (3.26) 

We note that A(t) is continuous for t bounded away from — oo. We need to show that 
our matrix A(t) is bounded as an operator on L l {t, oo) ® L 1 (t, oo) ® -^ 1 (t, oo) for this 
end we will use the Max norm. The entries of A(t) are ±q n (x) and ±p n (x). 

/ |g„(x)| dx = —= I \q{x) + /(x)n~^Gfe = Mi with M 1 < oo (3.27) 
it v2 i s 



We made use of the following change of variables together with formul4j(2.29) of [3] 

r— x i /— s , 

x = V2n + -p-r, and t = V2n + -p-r, (3.28) 
V2n6 V2n 6 

the fact that the asymptotics for p at infinity can be obtained from the following 
representation p = q + uq where it(oo) = or that u(x) is bounded for x away 
from minus infinity. We also assumed without lost of generalities for this section that 

2 3 

p(x) ~ Const • x 1 ' 4 e~3 x ^ as x — > oo which is a consequence of the asymptotics 

2 3 

q(x) ~ 2 ^ / _^ 1/4 e~3 x ^ as x — > oo. We also remarked that the scaled value of g„ in 
(2.29) of [3] is represented in terms of finite combinations of bounded functions (the 
itj's, Uj's, uVs ,) x\ i — 0, 1, 2, with p and g. So again we assumed that the scaled 

value of q n (x) was of order r nhx 1 e~^ x ^ 1 as x — > oo. A similar argument hold for p n , 
here we use formula (2.30) of [3] instead. 

/•oo /»oo 

/ = — / \q{x) + g(x)n~^ \dx = M 2 with M 2 < oo. (3.29) 

Jt v 2 J s 

Note that ||v4||Maz and ||J3||Afaa; are & t most 2Mi + 2M 2 . The fundamental local 
existence of solution for linear Ordinary Differential Equation says that equations 
have solutions in (a, oo) with a bounded away from infinity given by 



X(t) = exp f- jf" A(x)cfe^ • X(oo), and F(t) = exp (- J B(x)dx^j ■ Y(oo). 

(3.30) 

5 We set the constant c in (2.29) to zero, and use the known asymptotic of q(x) at infinity to 
deduce the existance of the integral. 
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This solution is convenient for the large n expansion of the probability distribution 
since it allows us to give a series expansion of the solutions of our solution in terms 
of q n and p n . The other advantage is the built in symmetries in matrices A and B. 
These symmetries make the computation of the matrix exponential very easy. We 
will start with the first system X{t) = exp (— J t °° A(x)dx) ■ X(oo). We set 



r / 



exp 



A{x) dx ] = exp < 



f™q n (x) dx \ ' 

— p n (x) dx 

A ir Pn(x)dx -f t co q n (x)dx 



= exp(M). 



M is of the form 



M 




/ J 

(3.31) 

(3.32) 



with 



2 k ~ 1 a k b k 



exp(M) 



Efc>i 



/ 1 + Efc>l (2k) 

_ 2 fc-l a fc-l fe fc + l 



2 fc-i a fe+i b fe-i 



V E 



(2fc)! 

2 fc ffl fc fe fc+i 
fe>0 (2fc+l)! 



E 



(2fc)! 

2 k - 1 a k b k 
(2k)\ 

2 k a k+1 b k 



E 



E 



2 fc a fc+i fe fc 

fc>0 (2fc+l) 
2 fc a fc 6 fc 



D! \ 



k>0 (2k+l)\ 



k>0 (2k+l)\ 

i _|_ V- 2 k a k b k 
1 l^k>l (2k)\ I 



We have 



and 



exp(M) n exp(M) 12 exp(M) 13 
exp(M) 21 exp(M) 22 exp(M) 23 
exp (M) 3i exp (M) 32 exp(M) 33 

u n ,e(t) = exp (M) 12 + c v exp(M)i 3 , 
V n>e {t) = exp(M) 22 + c^exp(M) 23 , 



Qn,e{t) = exp (M) 3 2 + exp(M) 



33 ■ 



(3.33) 

(3.34) 
(3.35) 

(3.36) 



3.2.1 Scaling 



At this point we scale the functions involved in (I3.33|) in terms of n at the point 
corresponding to the expected value of the largest eigenvalue. If we set 



t = r{s) = ^/2(n + c) + ^^, 

22?7,6 



(3.37) 



then equations (2.29) and (2.30) of [3] are 

Qn(r{s)) = Q n (r(s); r(s)) = ( q(s) + 



2c -1 



p(s) — cq(s)u(s) 



R3 
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(10c 2 - 10c + -)qi(s) + p 2 (s) + (-30c 2 + 10c + -)q(s)v(s) 
+p 1 (s)v(s) + p(s)vi(s) - q 2 (s)u(s) - q\{s)u x {s) - q(s)u 2 (s) 



+ (-10c 2 + -)p(s)u(s) + 20c 2 q(s)u 2 (s) 



n s 



+ 0(n~ 1 )e q {s) 



(3.38) 



and 



Pn{r{s)) = P n (r(s);r(s)) = n« q(s) + 



2c + 1 



p(s) — cq(s)u(s) 



(10c 2 + 10c + -)gi(s) +p 2 (s) + (-30c 2 - 10c + -)g(s)v(s) 

— _ 

+Pi(s)v(s) +p(s)v 1 (s) - q 2 (s)u(s) - gi(s)«i(s) - g(s)u 2 (s) 
3. 



+ -10c z 



■)p(s)u(s) + 20c 2 q(s)u 2 (s) 



_ 2 



20 



0(n- l )e p (s) 



If we change the variable in a and b by setting^] x := r(x), we obtain 



(3.39) 



v/2 



2c- 1 



p(x) — cq(x)u(x) 



77,3 



3 3 
(10c 2 - 10c + -)qi(x) + p 2 (x) + (-30c 2 + 10c + -)q(x)v(x) 

-pi(x)v(x) + p(x)vi(x) — q 2 (x)u(x) — qi(x)ui(x) — q(x)u 2 (x) 



+ (-10c 2 + -)p(x)u{x) + 20c 2 q{x)u 2 (x) 



and 



_ 2 
n 3 

~20~ 



+ O^^eJa;) dec, (3.40) 



a (s) + ax(s)n 1/3 + a 2 (s)n 2/3 + a 3 (s)n 1 



v/2 



q(x) + 



2c + 1 



p(x) — cq(x)u(x) 



77,3 



+ 



3 3 
(10c 2 + 10c + -)qi(x) + p 2 (x) + (-30c 2 - 10c + -)q(x)v(x) 

Pi(x)v(x) + p(x)v\(x) — q 2 (x)u(x) — qi(x)u\(x) — q(x)u 2 (x) 



■10c 2 + -)p(x)u(x) + 20c 2 q(x)u 2 (x) 



_ 2 
n 3 

~20~ 



+ 0(n~ 1 )e p (x) dx (3.41) 



a (s) + b^n- 1 ^ + b 2 (s)n- 2/3 + b 3 (s)n-\ 



6 We use the same letter in both sides in the change here to simplify notation. 
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We next focus on the following expression 

a k b k = (ab) k = (ajj + a (ai + &i)ra~s + (a (6 2 + a 2 ) + a^n'^ + Dn~ 1 ) fe . 
An expansion of this expression is 

a k b k = (ab) k = ajf + ka 2 ^' 2 (a (ai + bi)n~^ + (a (6 2 + a 2 ) + a^n'^ + Dn' 1 ^ + 

fc(fc-l) 



/ 1 2 

-aQ fc " 4 f a (ai + bi)n~* + (a (^2 + a 2 ) + ai&i)n~3 + Dn' 1 



2fc-2i 



i=3 



^a (ai + bi)n~? + (a (6 2 + a 2 ) + ai&i)n~s + Dn^ 1 ^ 
If we note that for i > 3 

(a (ai + + (a (b 2 + a 2 ) + ai&i)n~5 + Dn^ 1 ^ = 0(n _1 ), 

then the sum in this last term can be represented as 



+ af + kaf 2 + 



fc(fc-l) 2k- 



af~ 4 Oin- 1 ) =((al + l) k - af - kaf~ 2 



k(k-l) 
2 

k(k-l) 



,2fc-4 _ 2k _ hn 2k-2 



We have 



ka 



2k-2 




(2fc)! 



(a {ai + bi)n~5 + (a (6 2 + a 2 ) + ai&i)n~£ + DrT 1 j = (3.42) 



,2fc-l 
i 



2(2fe — 1)! 
and for k > 2 



(ai + &i)n 3 + (a 2 + 5 2 ) n s j + 



,2fc-2 
J 



2(2fc — 1)! 



a\b\n 3 



2k-2 
i 



2(2fc — 1)! 



k< 2(2k)\ ^ 4 ( a o( a i + ^i) n ° + ( a o( & 2 + ^2) + ai6i)n 3 + Dn 1 



Dn-\ 



(3.43) 



7 2fc-2 
*0 



7 2fc-2 
l 



S(2k-2)\ 8(2fc-l)! 
We have at this stage, 



(a 1 + b 1 ) 2 n--s + 0(n- 1 )y 



Ok-l„kuk ofc-l„2fc Qfe-l„2fe-l 

i+E ^ = 1 +E ? ?™3 L+ E 



fe>i 



(2fc)! ^ 2(2fc- 1)! 
fc>i v 7 fe>i v 1 



E 

fe>i 



r>fc-l 2fc-l 
z "0 

2(2fc - 1)! 



(a 2 + 6 2 ) 



+E 

fc>i 



fc-1 2fc-2 



2 fe - i ag- 
2(2fc- 1)! 



fc-1 2fc-2 



2 B - 1 a; 







ofc-l„2fe-2 

8(2fc-2)! ^8(2Jfe-l)! 



E 



(ai + 61 



n 3 + 



15 



fc>i 



(2*)! 



2^ 



fc>i 



(2fc- 1)! 



(3.44) 



(a 2 + & 2 ) ^v^o) 2 ^ 1 
2v ^ ^ (2*-l)l 



aA ^ (v^ao) 2 ^ 1 , / ^ (V2a ) 2k ^ (v^ao) 2 ' 



i(2*)! ^8 



n 3 + 



k>l 

(ai + 6 



2^ 



(2fc)! 

— sinh(V2ao) + 



(3.45) 



i^sinh(v/2a o) + 



— 5 (l + cosh(v / 2a c 

sinh(v / 2a-o) + ^ ^cosh(v / 2a ) - sinh(\/2ao) J (ai + &i) 2 

cosh(2ao + \)0{n~ x ) = exp(M)n = exp(M) 22 
A similar argument gives 



n 3 + 

(3.46) 



-E 



2fc _i afc+ i 6fc _ 



fe>i 



(2fc)! 



^ (l - cosh(v / 2ao)) + a \ bl (l - cosh(v / 2a )) 

2 \ / 2fln 



ai + &i 
2^ 



sinh(v / 2ao) 



n 3 + 



6i(6i-ai) a 2 -& 2 , /^i( a i _fe i) (ai + &i) 



2a 2 



2a r 



2a 2 



+ 



2a J K 0J I 16a 16a 8a 4 y v o; 



+ cosh(v / 2a ) 0(n _1 ) = exp(M) 12 , 
if we interchange a and b, then 



(3.47) 



E 

fc>i 



2 k-i a k-i b k+i 



ai + 6i 
2v/2 

a 2 - b 2 
2a 



(2k)\ 



sin 



- (l - cosh(A/2a ) 



bi — a\ 

2dn 



(l - cosh(v / 2a )) 



n 3 -f 



Qi(Qi ~h) | &2 - Q2 | ( Qi(&i - Qi) (ai + 6i) 2 



2a 2 



2a r 



2a 2 



) cosh(^ ) + f - ^ - ^ - ^1±M) sinh(V2a o) 

/ \ 16a 16a 8a 4 / 

+ cosh(v / 2a ) 0(n _1 ) = exp(M) 2 i, 



n 3 



We also have 



^ /rtu ■ 1NI = — sinh(\/2ao) + 



fc>0 



(2fc + l)! 72 



1 ~[ 1 cosh(v / 2a,o) + Ql t= 1 sinh(v / 2ao) 



2^00 



(3.48) 



n 3 + 



16 



V 8a 2a 2 ^ 1 o; V 4v^ 



(ai + bi 



8V^c 



2V2a 2 



+ 



a 2 - b 2 



J sinh(v / 2ao) 



n 



"i + sinh(v / 2a )O(n" 1 ) = exp(M) 



2V2a / 

and if we interchange a and 6 in this last formula, 
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-exp(M) 32 , 
(3.49) 



E 

fc>0 



2 k a k b k+1 



^ 2k + l y ^= sinh(V2a ) + — ^ &1 cosh(v / 2a ) + 7^7|~ sinh(V2a ) n 3 + 

/(a 1 + ^_^ az + bA j - / ( fll + M 2 _ (ax + 
V 8a 2a 2 ^ 1 0; V 4v^ 8v^a§ 2^ 



&2 - a 2 



] sinh(v2ao) 



"i + sinh(v / 2a )O(n" 1 ) = exp(M) 



31 



The last term of the exponential matrix (I3.33P is 

E2 k a k b k . r , ai + &i , . /- . 
/o; - t = cosh(V2a ) H j=— sinh(v 2a ) ri s + 



k>l 



(2k)\ 



s/2 



4 



-exp(M) 23 . 
(3.50) 

cosh(v^2ao) + 



( + ^ - i^l) sinh(V2a 



V\/2ao \/2 



4^0 / 



n s 



+ cosh(v / 2a )O(ra~ 1 ) = exp(M) 



33 • 



(3.51) 



We note that a/2oo is exactly the quantity \x defined in [23] as the new variable when 
solving for the limiting system of equation as n goes to infinity. 

We then use equations (I3.35p . (I3.46p . (13.501) together with the numerical value of 
c v for n even given by (12.221) to have the following expansion of V njB and g n e . 



VnM s )) = ^(l-exp(- v / 2a ))- 



a\ — b\ sinh(v / 2a ) Qi + b x 



2\[2 V2a 



2^2 



exp(-v / 2a ) 



n 3 -f 



(ai + b 



^-4 



\/2(a 2 + b 2 ) 



V2a ' 



+ 



2^00 



cosh(v^2cto) 



(Qi + ^i) 2 (1 



16a 2 , 



2^00 4a 



- 62) _ V2(a 2 + 62' 
4^2a 4 



\/2a ' 

sinh(v / 2a ) | n~*+0(-) 



and 



1 



q n ,e( T ( s )) = —^exp(-V2a ) + 



8 1 v^a 2a 2 



ai + 61 ^ ai - 61 sinh(v/2a ) 



v^ao 



n 3 



cosh(v^2cto) + 



^2 - «2 a 2 + &2 
. 2^00 2 



V^ax + fri) 2 | 1 ^ | fll 6i , ^(ai + ^i) 2 v^ 2 



16a^ 



4a^ 



sinh(v^ao) I n~f +0(-). 

/ n 
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3.2.2 Second system of equations involving the calligraphic variables for 
GOE n 



— — — n 

The system involving the calligraphic variables 

Y(t) = exp f- jf" B(x)dx\ ■ Y(oo). 



(3.52) 



We set 



roc 



exp ^— J B(x) dx 



Ai is of the form 



/ 



= exp < 
















k V -f t °°Pn{x)dx 

= exp(M). 



■ f t °° Qn(x) dx 



- J7° q n {x)dx \ 

- f™p n (x) dx 





/ J 



M = 



-a 
0-6 
-6 -a 



(3.53) 



with 



/ 1 I V 1 2 k - 1 a k b k 

1 "r Z^fc>l (2fe)! Z^fc>l 



exp(A^) = 



2 fc-i a fc+i fc fc-i 

(2fe)! 



E 2 fc a fc+l fc fc V 
fc>0 (2fe+l)! 



2 fc-i a fc-i fe fc+i 1 ^ 2 k ~ 1 a k b k _ 2 fc a fc fe fc +l 
2^fc>l (2fe)! 1 Z^fc>l (2fe)! Z^fc>0 (2fe+l)! 



_ 2 fc a fc b fc+1 _ 2 fc a fc+1 6 fc 1 i V ^ 

\ Z^fc>0 (2fc+l)! Z^fc>0 (2fe+l)! 1 ^ Z^fc>l 



2 k a k b k i 
fe>l (2fc)! / 



13 
23 
33, 



We have 



and 



(exp(M)n exp(A^)i2 exp(A^) 
exp(7W) 2 i exp(A4) 2 2 exp(A^) 
exp (M) 31 exp (M) 32 exp(A^) 

Q n;1 (t) = 2c tp exp(M) n + exp(.M) 13 , 
V n> i(t) = 2c v exp(M) 2 i + exp (M) 23, 

ft n ,i(t) = 2c v exp(A^) 3 i + exp(.M) 33 . 



(3.54) 



(3.55) 
(3.56) 

(3.57) 



We note that exp(_M 2 i) = -exp(M) 2 i, exp(_M 2 3) = exp(M) 23 , exp^^ 
— exp(M) 3 i and exp(.M 33 ) = exp(M) 33 . The solutions (13.561) and (I3.57P follow 
directly from the large n expansion obtained in the last subsection. We therefore 
the following solutions for V n ,\ and 7l n ^ 



have 



V n ,i{r{s)) = -J=(exp(-v / 2a )-l) + 



1 1 — 1 exp(— y/2a ) H — — - cosh("\/2ao) 



a 



V2i 
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b\ — a\ sinh(v / 2' 



a ) 



4a n 



2 

a 2 - b 2 
V2a 



+ 



n 3 + 



\/2a\ \f2a\b\ a 2 — b 2 



2al 



2a 2 , 



2a 2 



4a 
+ a 2 



cosh(v / 2a ) + 
-& 2 ~ 



4v/2 



2 v / 2a§ 



sinhfv^ao) ) n~* + O(-) 

n 



(3.51 



and 



7^ n! i(r(s)) = exp(-v / 2' 



a J 



ai — &i sinh(A/2a ) a\ + &i 
\/2 A/2a V2 



exp( 



a J 



n 3 



(ai + 6i) 2 a 2 + & 2 



4 

sinh^ao) + 



exp(— V2a 
a 2 - b 2 



4:V2a 



cosh(v / 2i 



a 



(gi ~ 
4V2a 



sinh(v / 2 1 



2a 



sinh(y / 2' 



a J 



^- cosh(\/2 l 



a J 



a 

2a; 



^ sinh(-s/2ao) 



o 



a 



_ 2 

n 3 



(3.59) 



3.2.3 Calligraphic variables for GSE„ 



We note that the GSE n case is identical to the GOE n up to a sign change and the 
parity of n for the calligraphic variables. The large n expansion for u n e (t) and v nje (t) 
follows from the matrix exponential f 1 3 . 3 3 j) . The boundary conditions need to be 
change to M nj£ (oo) = and u nj£ (oo) = and g nj£ (oo) = as n is odd. Therefore in 
this case 



u n ,e{t) 

Vn,e(t) 



exp(M) 12 
exp(M) 22 



(3.60) 
(3.61) 



and 



q n>£ (t) = exp(M) 32 (3.62) 

The large n expansions of these quantities is given by (I3.47p . (I3.46P and (I3.49P re- 
spectively. 

The system of equations satisfied by the calligraphic variables is 

, / Qn A (t) \ ( -q n (t) \ ( Q nA {t) \ 

-r. V nA {t) = - Pn (t) ■ V nA (t) . (3.63) 

M V n nA (t) J \ - Vn {t) -q n (t) o / V K nA (t) ) 

where lZ nA (t) = 1 + 1Z nA (i). The boundary conditions in this case are 

Qn,4(°°) \ / -c M \ / 



V nA {oo) 
K nA (oo) 



-C0 



as n is odd 



(3.64) 



We can use the same technique as for the Orthogonal case to find a large n expansion 
of Qn,4(*), Vn^t) and U nA {t). 
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3.3 Large n expansion of the probability distribution of the 
largest eigenvalue 



We recall that the quantity of interest is (12.241) . Under our change of variables it 
reads 



(l-v nte )(l--n n>1 )-hq ntE -c tp )V n ,i = ^ 



V n;S (l + n n<1 ) - V n ,i(q n , £ - Cp) . (3.65) 



Upon substitutions of the newly derived expressions in the right of (13.651) . the right 
side of (I3.65P takes the form 



Ql n bl (l -exp(-v / 2a )) - ^ + J >1 exp(- v / 2a ) 



exp(-v / 2a ) + 

aibx - a\ t a 2 -b 2 , ( 3a 2 - b\ (ai + bi) 2 a 2 + b 2 

lV2a 16V2O0 4V2 



n s 



2a r 



+ 



exp(— 2v^2a-o) 



+ 



(ai + 6i) 2 3(a 2 + 6 2 ) «i («i + h 



+ 



+ 



Ay/2 2a 2 , 16v^ao 2v^ao 4v / 2a 2a o 
(«! + b] 



exp(-v / 2a ) 



— sinh(\/2ao) - cosh(v / 2a ) + O(-) (3.66) 
We follow Tracy and Widom [21] and denote by fi(s) the quantity \^2a (s), 

POO 

(j, := /jl(s) = V2a {s) = / q(x)dx, (3.67) 

J s 

and we introduce the following notations, 

) POO 

p(x)dx, a := a(s) = / q(x) u(x) dx, (3.68) 
1 r°/2c-l 



^ := v[s] 



2 

2c + 1 



p(x) — cq(x)u(x) J dx, 
p(x) — cq(x)u(x) ] dx, 



W ' ) = ^f ( 

If 00 / 3 3 

« 2 (s) = — ^ y f (10c 2 - 10c+ -) 0l (x) +p 2 {x) + (-30c 2 + 10c+ -)q( 



x)v{x) + pi^(x) 



-p(x)t>i(x) — q 2 (x)u(x) — q\(x)ux{x) — q{x)u 2 {x) + (- — 10c 2 )p(x)u(x) + 20c 2 g(a;)u 2 (x) J dx, 



If 00 / 3 3 

62(a) = J ^(lOc 2 + 10c + -) qi (x) + p 2 (x) + (-30c 2 - 10c + -)q{ 



x)v(x) + p\v{x) 



20 



3 \ 
-p(x)vi(x) — q 2 (x)u(x) — qi(x)ui(x) — q(x)u 2 (x) + (- — 10c 2 )p(x)u(x) + 20c 2 q(x)u 2 (x) J dx. 

We note that 

1 f°° 

ai(s) — bi(s) = ■= I p(x)dx, (3.69) 

v2 Js 

1 r°° 

ai(s) + bx(s) = —j= I (2cp(x) — 2cq(x)u(x))dx = — V2cq(s), (3.70) 
v 2 Js 

c f°° c 
a 2 (s) - b 2 (s) = -j= J (q(x)v(x) - q x {x))dx = ~^P( S )- ( 3 - 71 ) 



and 



We set 

20^2 



r)(s) = a 2 (s)+b 2 (s) = — ^ I ( (20c 2 +3)q 1 (x)+2p 2 (x)+(-60(?+3)q(x)v(x)+2p 1 v(x) 



+2p(x)vi(x) — 2q 2 (x)u(x) — 2q 1 (x)ui(x) — 2q(x)u 2 (x) + (3 — 20c 2 )p(x)u(x) 

+ 40c 2 g(x)u 2 (x)^ dx. (3.72) 

20c 2 q'(s) + 3p(s) 1 



+ 



(6qv+3pu+2p 2 +2p 1 v+2pvi— 2q 2 u— 2q\U\ — 2qu 2 ) (x) dx. 

(3.73) 



20V2 20^2 
This last equality comes from 

(20c 2 + 3)q l (x) + (3 - 60c 2 )q(x)v{x) + (3 - 20c 2 )p{x)u{x) + 4Qc 2 q{x)u 2 {x 



20c 2 [q\ (x) — 3q(x)v(x) — p(x)u(x) + 2q(x)u 2 (x)) + 3(qi(x) + q(x)v(x) + p{x)u{x)) , 

(3.74) 

the substitutions 

qi(x) = xq(x) — q(x)v(x) + p(x)u(x), and qi(x) = p'(x) + q(x)v(x), (3.75) 
in jXH]) 

20c 2 (xg(x) — 4q(x)v(x) + 2g(x)w 2 (a;)) + 3(p'(x) + 2g(x)u(x) 

= 20c 2 (xq(x) + 2q(x)(-2v(x) + u 2 (x))) + 3(p'(x) + 2q(x)v(x) +p(x)u(x)), (3.76) 
and the substitutions 

q 2 (x) = u 2 (x) — 2v(x), and q" (x) = xq(x) + 2q 3 (x) (3.77) 

in fl3~76|) 

20c 2 (xg(x) + 2g 3 (x)) + 3{p'(x) + 2q(x)v(x) + p(x)u(x)) 

= 20c 2 q"{x) + 3p'{x) + 6q(x)v(x) + 3p(x)u(x). (3.78) 
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With these representations, (13.661) is 



cg(s)e" M - — -e" M ) 



c 2 - 2c + 1/4 2 c(c- 1) 



n 3 -f 



2c- 1 2 c r° . . . cp(s) 
v — ——^v I q[x)u[x)dx + 



8/i 



-i/ — 



V 2/x 2 y s " v ' v ' 1 2/i 

9 / /-00 



i/ / q(x)u(x)dx H 

4/x y s 8/i 



q(x)u(x)dx 



2 „2 „2 



c 2 g 2 (s) 
8// 



4V2. 



-2/i 



s ^/-" 

r c 2 g 2 ( S ) 377 2 -// /c 2 -c+ 1/4 2 2c 2 - c r° 

1/ — i/ / q{x)u{x)dx 



2 4v/2 2/i 2 



,2 /-co 



+ "^"( / q{x)u{x)dx)' 



c 2 q 2 (s) c 2 -l/4 
8/i 8/t 



-— ( / q(x)u(x)dx 



c 2 g 2 (sj 



z/ 2 H i> / q(x)u(x)dx 

4/i 



c 2 - 1/4 



cp(s) 
2/i 



2/i 2 



e M H — ^ ^ 7 sinh(/x)— I ' z/ 2 — c 2 ^ / q(x)u(x)dx 



( / g(x)u(a;)(ix)' 



2^ cosh(/t) I 



n"3 + O(-) 



as n goes to infinity uniformly for s bounded away from minus infinity. 



(3.79) 



Finally we combine (13.791) with Theorer dl.il to have the following version of The- 
oren il.2l 



If we set t = r(s), then as n — > oo 
F 2 (t) = F 2 (s)-L-» + 



c(q(s) + u(s))e-»-—(l-e-») 



1 



ca(s) cp(s) (2c — 1) v 2 



4/i 2 



n s + 



C «(a)(cg( a ) e-" _ (i _ e -/»; 



+ e 



- 2 , / , c 2 « 2 ( g ) _ c 2 g 2 ( g ) _ (c 2 -c)^(s) (| - 2c + c 2 ) 



4^ 8/x 



8/t 



4/t 



8/i 



^fc 2 q 2 (s) 3t] c 2 a 2 (s) cp(s) c 2 g 2 (s) c 2 v a(s) {-\ + c 2 )v 2 



Ay/2 4/i 2/t 8/i 



4/i 



8/i 



2-/i / c 2 a 2 (s) _ (2c 2 -c)va(s) (\-c+ c 2 ) v 2 
2/i 2 1 



c 2 a 2 (s) — c 2 v a(s) 



(-J + c 2 ) z/ 2 \ cosh(/i) c 2 g 2 (s) 



+ 



8/i 2 



sinh(/i) 



n 3 j> + 0(n x ) 
(3.80) 



uniformly in s. 



To simplify the n~3 term in (I3.80p . we use the representation p(s) = q'(s)+q(s)u(s) 
which says in this setting that u(s) = a(s) — q(s). The result of this substitution is 
Theorem ! 1.21 
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4 Conclusion 



We note that unlike F n ^{t) for the GUE n , the GOE n large n expansion of the proba- 
bility distribution of the largest eigenvalue F n> i(t) has a non vanishing correction 
term. Thus the convergence to the limiting Tracy- Widom distribution F\(t) is slower. 
Numerical applications of F n ^{t) follows easily from q(s) this is one consequence of 
our representation of F n! i(t) in Theorer dl.21 All the terms on the right side of (11.361) 
can be expressed in terms of q(s) and q'(s). 

The GSE n largest eigenvalue distribution is derived in a similar way (the only major 
difference being that n needs to be odd in this case.) 
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